Fractionation of polydisperse systems: multi-phase coexistence 
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The width of the distribution of species in a polydisperse system is employed in a small- variable 
expansion, to obtain a well-controlled and compact scheme by which to calculate phase equilibria in 
multi-phase systems. General and universal relations are derived, which determine the partitioning 
of the fluid components among the phases. The analysis applies to mixtures of arbitrarily many 
slightly-polydisperse components. An explicit solution is approximated for hard spheres. 

PACS numbers: 05.20.-y, 64.10.+h 



It is vital to gain an understanding of polydispersity, 
due to its ubiquity in both synthetic and biological com- 
plex fluids. A polydisperse substance is a mixture of in- 
finitely many components, and can, in general, separate 
into arbitrary numbers of coexisting phases. These prop- 
erties typically engender great mathematical complexity, 
and have been a stumbling block to the concise formu- 
lation of polydisperse thermodynamics. Experimental 
and simulational QJ^] studies of polydisperse poly- 
meric fluids and colloidal suspensions have catalogued 
diverse behavior and intricate phase diagrams. Until re- 
cently, theoretical treatments of polydispersity relied on 
uncontrolled approximations M or idealised models , 
while generic schemes and fundamental understanding 
remained elusive. 

The phase behavior of pure (i.e. monodisperse) systems 
is (in principle at least) relatively straightforward to anal- 
yse. The standard method, formulated last century ||, 
involves integrating (by various approximate methods) 
the Boltzmann factor over all configurations to construct 
the Helmholtz free energy as a function of temperature, 
density and volume. From this, the densities of coexisting 
phases can be calculated and the phase diagram deduced. 
One source of difficulty in analysing the phase equilibria 
of polydisperse systems is that the density alone docs 
not fully characterize a phase. Instead, we wish to know 
the entire composition of each phase. That requires the 
evaluation of an infinite set of variables. 

Two systematic schemes were developed recently to 
solve the polydisperse phase equilibria problem. The 
powerful 'annealed moments' method of Sollich and 
Cates |l(J and Warren H applies to a large subset of 
model systems, and will not be further discussed here. 
The second scheme, which applies to real systems, was 
developed by the present author |fl2| . It uses the width 
of the distribution of species as a small expansion pa- 
rameter, and is therefore valid for slightly polydisperse 
systems. In other words, the scheme is applicable when- 
ever the polydisperse property (e.g. the size or charge of a 
particle) varies only a little throughout the system. (N.B. 
Here, 'particle' is used to denote the polydisperse fluid 
elements, which could be polymer molecules, colloidal 
latices, etc.) The method was used [Q to find the com- 
plete distributions of species in two coexisting phases of 



any slightly polydisperse system, resulting in a universal 
law of fractionation. In this paper, the method is applied 
to coexistence between arbitrary numbers of phases — a 
situation of importance to many polydisperse substances 

II- 

A slightly polydisperse system (one with a narrow dis- 
tribution of species) is in principle very different from a 
truly pure one, which has no mixing entropy and whose 
distribution is a Dirac delta function. Nevertheless, one 
would expect the physical properties of the two systems 
to be very similar. That similarity motivates this study, 
since the pure system is vastly simpler to analyse than 
its polydisperse counterpart. To exploit that simplicity, 
a formalism is required which treats mono- and poly- 
disperse systems on an equal footing. Such a formalism 
is now derived (extending the method for two-phase equi- 
librium in Ref. Jl|). 

Let us first define a number Si to characterize each of 
the N particles in the system (with i = l,...,N) [p"4[ . 
For size polydispersity, this is the fractional difference 
£j = (Ri — R )/R of the particle's radius Ri from some 
reference length Rq (with the obvious generalization to 
charge polydispersity etc.). Henceforth, e shall be re- 
ferred to as the size parameter, for definiteness. The 
population of species in the system is characterized by a 
continuous distribution f(s), which is unnormalized so 



f(s)de = N. 



In general, the free energy F of a polydisperse system 
is a complicated functional of f(e). It will be expressed 
in units of UbT where ks is Boltzmann's constant and 
T is temperature. For a polydisperse ideal gas, the free 
energy F ld is easily shown [|7|Jll|] to be, per unit volume, 
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which is the usual ideal gas free energy, summed over all 
species. As this expression contains the mixing entropy, 
it is useful to write the free energy of a non-ideal system 
as 



F = F id 



(2) 
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Here F ex is the 'excess' part of the free energy (over and 
above the ideal part), deriving from interactions. 

Let us consider a system whose 'initial' population (be- 
fore phase separation) is known. This will be called the 
'parent' distribution fp{s). In a system where this parent 
is partitioned into Ai coexisting phases, we wish to deter- 
mine the distribution f(e) A in each phase A. = 1, ...,M.. 
By conservation of matter, 
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A=l 



At equilibrium, the chemical potential is equal in all co- 
existing phases. This statement applies for each species 
of particles, so the equation 



P&A = Me)f 
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represents an uncountable infinity of thermodynamic 
constraints for any pair of phases A and B. Since there is 
a continuum of species, the chemical potential is a func- 
tional derivative of the free energy 



p,{e) = 
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From Eq. |^, fi(e) can be written in two parts 

M ( £ )=/i id ( £ )+/x ex ( £ ). 
Functional differentiation of Eq. |l| yields 



fi id (e) =ln 
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Collecting together Eqs. ||, || and ^ gives the ratios of 
densities in any two of the A4 coexisting phases 
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in terms of the excess parts of the chemical potentials. 
Thus, all but one distribution can be eliminated from 
Eq. ||, yielding the solution for any given phase 



f(e) A = f P (e) 



B=l 



exp {^{e) A -^(e) B ) (9a) 



where A* ex (e) 
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Given a knowledge of F cx , which specifies the interac- 
tions in the system, (and of the phase volumes), Eqs. |9a| 
and ^b| represent a complete solution to the problem. 
However, they constitute an uncountable infinity of non- 
linear simultaneous equations. This is the source of the 
mathematical complexity mentioned earlier. 



Some simplification is achieved by making a change 
of variables. Rather than expressing a thermodynamic 
state in terms of the densities of the individual species of 
particles, f(s)/V , let us use moments of this distribution 
(as in [noUOJlO]). The thermodynamic variables 
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a = 0, 1, oo 
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will be called 'moment densities'. Note that p a = e a p, 
so that po is the overall number density p. [Mean powers 
of the size parameter, e a , are moments of the normalized 
distribution p(e) = f(e)/N.] Each moment density, be- 
ing a linear combination of conserved species densities, is 
itself conserved and, accordingly, respects the usual equi- 
librium conditions. For instance, each 'moment chemical 
potential', defined by p a = d(F/V)/dp a , is equal in co- 
existing phases. This is clear from expanding the species 
chemical potential in partial derivatives 
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Thus, equality of p(e) in coexisting phases requires equal- 
ity of each p a . 

We now have a discrete set of thermodynamic vari- 
ables, and can substitute the power series expression 
(Eq. |ll|) for p(e) into Eq. |a], yielding 



/ B=l \a=0 



(p c a A-PT B )e a 



with 
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which, with Eq. |l^, form a countable infinity of simulta- 
neous equations. The excess free energy is now & function 
F ex (po, pi, ...) of the moment densities. 

The equations thus far are perfectly general, but the 
advantages of this formalism become apparent when we 
consider a narrow distribution of sizes, i.e. a system 
which is close to monodisperse. If the origin for the pa- 
rameter e is chosen (by fixing the reference i?o) 1° be 
close to the centre of the narrow distribution, then e is 
a small number for most if not all particles. Hence in 
Eq. 12a, fp(e) vanishes for large e, so the power series 
in the denominator becomes a well-controlled expansion. 
The results have a particularly simple form if the origin 
is chosen to be the mean of the parent distribution, so 
that ep = 0. Henceforth t his c hoic e is a ssumed. 

The solution to Eqs. ^0. 12a and 12b is now calculated 
to first order in e. This will yield the exact phase equilib- 
ria i n th e limit of a narrow parent, e 2 p — > 0. Expanding 
Eq. 12a to first order and integrating over e gives 
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To zeroth order, Eq 
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Note the different ord ers o f expansion. Substituting these 
expressions back into 12a yields 



f(e) A = fp(e) 
N A N 
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To obtain Eq. |13|, the pref actor of fp(e) in Eq. |l2a| is 
expanded to first order in e, but the distribution fp(s) 
itself remains exact. Thus, other than narrowness, no 
limitations are put on the form of fp(s). Any distribu- 
tion can be treated, however asymmetric or discontinu- 
ous. The population may even contain finite amounts of 
some components, contributing delta functions to fp(e). 

In Eq. 13 we see that the distribution in any given 
phase A depends, as one would expect, on the proper- 
ties of all the other M. phases with which it coexists. 
However, taking the difference (denoted A) between the 
normalized distributions in any two of the M coexisting 
phases, we find the strikingly simple expression 



Ap(e) -> -epJe) 



(14) 



in the limit as e 2 p — > 0, where p p (e) is the normal- 
ized parent distribution. [Note that the solution for each 
phase (Eq. |l^) is recoverable from the neater sum (Eq. |^) 
and difference (Eq. |l4|) equations.] Surprisingly we have 
found that, in the multi- phase system, the difference in 
compositions of any pair of phases is identical to the 
expression found earlier [l2| for two-phase coexistence. 
Thus the same universal laws follow pl| , relating any 
pair of phases. This is not an obvious result, since the 
parent appearing in Eq. [m| is the combined population of 
the whole system, not just of the two phases in question 
as it is in the two-phase coexistence problem. 

Equation [ti] is very generally applicable. It is valid for 
any system with a narrow distribution (that is, narrower 
than the range of lincarizability of the fugacity), what- 
ever particles or interactions it comprises. Furthermore, 
recall that e need not parameterize size deviations, but 
could represent charge, mass or any other sole polydis- 
perse quantity. By analysing multi-phase coexistence, we 
have found that Eq. |l4| does not even depend on M. , the 
number of phases present. 

We have considered a system in which a slightly poly- 
disperse fluid component is partitioned among several 
phases. The coexistence of more than two phases may be 
the result of tuning the temperature to the triple point 
of the monodisperse reference system. Alternatively, the 
slightly polydisperse particles may be in the presence of 
other, dissimilar components which, by the Gibbs phase 
rule (which states that an n-component mixture can ex- 
hibit up to n + 1 coexisting phases at arbitrary tempera- 
ture), can induce multi-phase coexistence [[17|. Within 



such a multi-component system, a particular, slightly 
polydisperse component will respect the above relations, 
which may be tested by an experimental probe which 
is 'blind' to the other components. For instance, near- 
monodisperse colloidal particles in the presence of 'de- 
pletant' species jl8| exhibit multiple phases. Light scat- 
tered only from the near-monodisperse colloid contains 
information on its fractionation [p.9| , which should obey 
Eq. [l4| As an illustration, a multi-phase colloidal sample 
with the composition shown in the figure, will obey the 
above relations, applied only to those particles in range 
X, with the origin of e defined at its centre. The rela- 
tions are equally applicable to particles in range Y, if we 
are blind to all other particles (e.g. they may be made 
invisible by matching their refractive index to that of the 
solvent), and redefine e = appropriately. 

The form of the solution in Eq. [li] is of interest in it- 
self, not least for the non-appearance of A4. However, 
one quantity remains unknown: the constant of propor- 
tionality A/i^ x . That constant is system-dependent. For 
some substances, /if x can be calculated using thermody- 
namic perturbation theory [jl2 15 1. Unfortunately, this is 
not possible for a system of hard-spheres, as it's Hamilto- 
nian is non-differentiable. Since the hard-sphere system 
is of great practical interest for modelling systems with 
repulsive interactions, the constant of proportionality is 
now calculated for that case. 

The excess part of the free energy of the polydisperse 
hard-sphere fluid can be Taylor expanded in the small 
size parameter of each of the N particles of interest thus 
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where F™ no 



is the excess free energy of the reference 
component of monodisperse hard spheres (in the pres- 
ence of the rest of the system — see Fig. [[]). In the 
reference component, all particles are alike, so the differ- 
entiation may be performed on particle number 1 only, 
without loss of generality, giving 
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£1=0 



The change in the identity (the species) of particle 1 when 
its size is varied affects only F ld . The excess free energy 
contains the physical effect of the particle's size on the 
rest of the system. By its presence in the container, parti- 
cle 1 simply excludes other particles from a volume V exc i, 
given that its interactions are purely hard and repulsive. 
Thus, increasing its size reduces the effective system vol- 
ume, so 
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In fact the volume from which particle 1 excludes other 
particles, V eKC \, depends on their species, so the quantity 
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in Eq. [15] is a net effective value, defined by the equa- 
tion. For the special case of an almost pure hard sphere 
system (not in the presence of other, dissimilar compo- 
nents), Voxel = §7!\ff P (2 + £i) 3 at low density (correct 
up to second virial coefficient). At high density, the ge- 
ometry of high-order inter-particle interactions modifies 
this. In any case, dV cyic \/ de\ is of order a particle volume. 
The resulting excess free energy density of a polydisperse 
hard sphere fluid is 
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where V oS is some (unknown) effective volume, of order 
the volume of an average sphere, and exactl y tha t for a 
near-pure, low-density system. Applying Eq. 12b yields 



Hf = 12P ex V e 



(17) 



in terms of the system's excess pressure P ex over an ideal 
gas. Since coexisting phases have the same total pressure, 
it follows that A/xf x = -12 V cS AP id . So the constant of 
proportionality in Eq. [li] is 



T2V ott Ap 



(18) 



for hard spheres in ergodic (fluid) phases. This calcu- 
lation contains the lowest-order effects of polydispersity. 
Once the polydispersity is sufficient to alter the mode of 
packing (e.g. small particles preferentially filling the gaps 
between big ones), higher-order analysis is needed. 

It is apparent that combining a moment descrip- 
tion with a small-variable expansion in the distribu- 
tion's width is a productive way to analyse polydis- 
perse systems. While the applications of this study 
are clearly wide-ranging, it is intended to extend its 
scope by analysing correlation functions and multiply- 
polydisperse systems fl5| . In addition, some work is re- 
quired, using higher-order analysis, to establish the ra- 
dius of convergence of the expansion, and quantify more 
precisely the method's regime of validity. 
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FIG. 1. The composition of a sample which could exhibit 
multi-phase coexistence. The narrow part of the distribution 
in range X can be treated by the present theory, as could the 
part in range Y. 
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